Journal of Nonlinear Mathematical Physics Volume 8, Number 2 (2001), 256-287 Article 



o 



> 



Reflectionless Analytic Difference Operators II. 
Relations to Soliton Systems 



S N M RUIJSENAARS 

Centre for Mathematics and Computer Science 
P.O. Box 94079, 1090 CB Amsterdam, The Netherlands 

O \ Received July 30, 2000; Accepted January 10, 2001 

(N 
< 



Abstract 



This is the second part of a series of papers dealing with an extensive class of ana- 
lytic difference operators admitting reflectionless eigenfunctions. In the first part, the 
pertinent difference operators and their reflectionless eigenfunctions are constructed 
from given "spectral data" , in analogy with the 1ST for reflectionless Schrodinger and 
Jacobi operators. In the present paper, we introduce a suitable time dependence in 
^0 \ the data, arriving at explicit solutions to a nonlocal evolution equation of Toda type, 

pH . which may be viewed as an analog of the KdV and Toda lattice equations for the 

latter operators. As a corollary, we reobtain various known results concerning reflec- 
tionless Schrodinger and Jacobi operators. Exploiting a reparametrization in terms 
of relativistic Calogero-Moser systems, we also present a detailed study of A^-soliton 
solutions to our nonlocal evolution equation. 

> 

m 

■ 1 Introduction 

O 

\ In our previous paper Ref. we have introduced and studied an extensive class of AAOs 

(analytic difference operators) admitting reflectionless eigenfunctions. Our construction in 
Ref. is patterned after the 1ST for Schrodinger and Jacobi operators. In brief, we start 
^ \ from given spectral data {r,fj.), and associate to these data a meromorphic reflectionless 

wave function W{x,p) and meromorphic coefficients Va{x), Vb{x) ("potentials") of an AAO 
A of the form 



A = Ti + Va{x)T^i + Vb{x), T±i = expi^id,), (1.1) 

such that we have the eigenvalue equation 

{AW){x,p) = {f + e"P) W(x,p). (1.2) 

The data r = (ri, . . . , r^r) allowed in Ref. Q (from now on denoted by Part I) gives rise 
to the poles of the transmission coefficient, just as for Schrodinger and Jacobi operators. 
It consists of complex numbers restricted by 

Imr„ E (-7r,0)U(0,7r), n = l,...,iV, (1.3) 
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and 

gr„_^g±r„^ 1 < < n < iV. (1.4) 

This guarantees that the Cauchy matrix 

Cir)mn = — — - — m,n = l,...,N, (1.5) 

g' m — n 

is well defined and regular, just as the Cauchy matrix pertinent to the Schrodinger and 
Jacobi cases. (Cf. Eqs. (2.30)~(2.33) in Part I, or, briefly, I(2.30)-(2.33).) 

The "normalization coefficients" fj, = (^i, . . . , hn) permitted in Part I are however far 
more general than for reflectionless Schrodinger and Jacobi operators. Indeed, they are 
allowed to be meromorphic functions satisfying (cf. 1(2.34)) 

Unix + i) = fJ-ix), lim /i„(x) = c„, c„ G C*, n=l,...,N. (1.6) 

|Re >oo 

Of course, this includes the constant case jU„(x) = c„, which is the analog of the 
Schrodinger and Jacobi settings. In Part III of this series of papers (which deals with 
various functional-analytic features), we focus attention on the constant multiplier case. 
In Sections 2 and 5 of the present paper, however, we allow the same multipliers as in 
Part I. 

In Section 2 we introduce time dependence in the multipliers fii{x), . . . , fi]y{x). Ac- 
cordingly, the potentials Va{x), Vb{x) and wave function W{x,p) depend on time. The 
time dependence is chosen such that for each set of "initial" spectral data (r, fi), we arrive 
at a solution to a nonlocal Toda type equation. Specifically, this equation reads 

^{x, t) = i exp(i[^'(x + i,t) - ^'(x, t)]) - i exp(i[^'(x, t) -^{x -i, t)]). (1.7) 

Here, the function exp[i^(x, t)], x, t E C, is assumed to be meromorphic in x for all t. With 
suitable restrictions on the data r, /u, we show that our solutions to (|1.7| ) are real-valued, 
real-analytic functions for real x, t, with a solitonic long-time behavior (in Section 6). 

Soliton equations of this nonlocal type have been introduced and studied before, cf. San- 
tini's review Ref. and references given there. (In particular, when one replaces the Ihs 
of (0) by idx^{x,t), then one obtains the so-called intermediate Toda lattice [^.) Even 



so, it seems that (1.7) is a new soliton equation, in the sense that there appears to be no 
obvious transformation relating it to previously known soliton equations. 

We should stress at this point that we do not associate to (|l]^) a clear-cut auxiliary 
linear spectral problem for self-adjoint AAOs on L^(R, dx). In fact, we view it as a 
challenging open problem to do so, in such a way that the self-adjoint AAOs of Part III 
arise as the ones with vanishing reflection. 

In Sections 3 and 4 we present further results that render the existence of such a scenario 
plausible. Indeed, in these sections we show that our (constant multiplier) AAOs can be 
related to the well-known reflectionless Schrodinger and Jacobi operators, respectively. 



Moreover, in Section 4 we show that the corresponding solutions to (1.7) obtained in 
Section 2 are related in the same way to the Toda lattice solitons. 

Just as Part I Q, Sections 2-4 are basically self-contained, in the sense that we need not 
invoke any substantial previous results in the literature. (To be sure, the previous, "IST- 
based" , literature in the KdV and Toda lattice settings did provide considerable inspiration 



258 



S N M Ruijsenaars 



for this paper as well as Part I.) For example, we obtain the solution property for the Toda 
lattice solitons as a corollary of the corresponding result in Section 2. 

In contrast, Section 5 (and Section 6, too) involves previous work of ours: We demon- 
strate there that the reflectionless wave functions and AAOs can be connected to a special 
class of relativistic Calogero-Moser type systems. We have labeled the pertinent finite- 
dimensional soliton systems by IIrei(7" = l'^ in our paper |^]; more general integrable 
A^-particle systems of Calogero-Moser and Toda type are surveyed in our lecture notes 
Ref. §. 

This relation plays an important role in Part III. It enables us to invoke various results 
from Ref. to control analytical difficulties. Staying within the context of the present 
paper, it yields a map from an arbitrary point in the Ilre^T" = '''"/2) A^-body phase space 
to a (real- valued) A-soliton solution of ( |1.7D . As such, it gives rise to one more example 
of what we have dubbed soliton-particle correspondence in our survey Ref. In this 
connection we also mention previous results on this correspondence for the KdV and Toda 
lattice solitons [^], and more recent results on the relation between relativistic Calogero- 
Moser systems and the 2D Toda field theory |^ IC]. 

A novel feature of the correspondence in the present setting is that our general A-soliton 
solutions to (|1.7] ) are encoded via the defining Lax matrix of the IlreKT" = '^1'^) system, 
and not via the dual ("action- angle") Lax matrix (as is the case for the sine-Gordon 
and modified KdV particle- like solutions |5|). To be more specific, the general A-soliton 
solution consists of A+ G {0, 1, ... , A} solitons moving to the right, and A_ = N — 
solitons moving to the left. The right-movers and left-movers are parametrized by the 
particle and antiparticle variables, respectively. When A_|_ or A^„ vanishes, one is dealing 
with a self-dual pair of Lax matrices, so this new feature plays no role. 

For A_(_A_ > 0, however, self-duality breaks down. In that case the correspondence 
between right-/left-movers and particles/antiparticles is quite different from the corre- 
spondence between solitons, antisolitons and breathers in the sine-Gordon and modified 
KdV settings on the one hand, and particles, antiparticles, and their bound states in the 
IIi-ei(T = 7r/2) system on the other hand. 

The consequences of this novel type of correspondence are made explicit in Section 6, 
where we study the A^-soliton solutions to (1.7). First of all, we demonstrate that the 
pertinent solutions deserve their name. Indeed, we show that for long times they can be 
approximated by linear combinations of A 1-soliton solutions; moreover, the asymptotic 
velocities are conserved and the position shifts are factorized in terms of pair shifts. 

We obtain the features just mentioned in a quite direct and elementary fashion, yielding 
however a weaker approximation result than we obtained for the sine-Gordon and (m)KdV 
solitons in Section 7 of Ref. Just as in all cases studied previously, the long-time 
behavior of the A-soliton solutions is intimately related to the spectral asymptotics for 
t — > zboo of time-dependent matrices. When A_ or A+ vanishes, we show that our results 
in loc. cit. give rise to a natural notion of global, non-intersecting, space-time trajectories 
for the A^ right-moving or A left-moving solitons. 

For A+A^_ > 0, however, we wind up with non- intersecting space-time trajectories only 
for sufficiently large times. In between, trajectories are ill-defined, since a collision of two 
trajectories typically gives rise to a complex-conjugate pair of eigenvalues. From a physical 
viewpoint, the right-moving soliton ("particle") and left- moving soliton ("antiparticle") 
involved in the collision form a resonance for a certain period of time. 
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2 A nonlocal Toda type soliton equation 

We begin by collecting the definitions of various important quantities, cf. Section 2 in 
Part I. We have already recalled the restrictions on the spectral data (r, fi) and the defi- 
nition of the Cauchy matrix C(r), cf. (pr^)-( |L6| ). The dependence on is encoded in the 
diagonal matrix 

D{r, fi; x) = diag {d{ri, fii; x), d{rN,fJ'N] x)), (2.1 
where the function d is defined by 

j u{x)e-^^P-, ImpG(0,7r), 

To ease the notation, we often write 

(i„(x) = d(r„,/i„;x), n = l,...,N. (2.3 

All of the remaining quantities can now be defined in terms of the solution R{r, /i; x 
to the linear system 

{D{r,fi;x) + C{r))R{x) = C, C = (1, • • • , !)*• (2-4 
Specifically, introducing the auxiliary functions 

N 

X{r,fi;x) = l + ^e'^'^Rn{r,fi;x), (2.5 

n=l 

N 

S(r, //; x) = ^ Rn{r, /i; x), (2.6 
n=l 

the potentials are given by 

Va{r,fj,;x) = X{r,fi;x)/X{r,iJ,;x (2.7 

Vb{r, fi; x) = T.{r, fi;x -i) - S(r, fi; x), (2.8 

and the wave function reads 

^ir,,;.,p)^e-^(l-f:^^). (2.9 

\ n=l / 

We now introduce time-dependent multipliers 

//„(r„;x,t) =^n(x)exp(ii [e''" -e"''"]), n = l,...,iV. (2.10 



Correspondingly, the above quantities (2.1)~( ^^ are henceforth viewed as depending on 
time as well. But as a rule we suppress the dependence on t, just as the dependence on 
(r, /i). Until further notice, t may be viewed as a complex parameter. We denote partial 
differentiation with respect to t by a dot. (The factor i in ( 2.10| ) occurs with an eye on 
later reality restrictions.) 

We proceed by obtaining the time derivatives of the above quantities. The pertinent 
formulas are far from immediate, so for clarity we assemble them in a series of propositions. 
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Proposition 2.1. One has 

R{x) = i diag (e'^'i , . . . , e"'^'^) {R{x) - R{x - i)) - iVb{x)R{x). 
Proof. Combining (|2l|), (U) witli (|2.10D , we obtain 

D{x) = iD^D{x), 
wliere we have introduced 

D_ = diag {e''' - e-''^ ,...,6'^ - e"''^) . 
Now from ( |2.4D we deduce 

D{x)R{x) + {D{x) + C7)i?(x) = 0, 
so that 

(Z)(x) + C)R{x) = -iD^D{x)R{x). 
Therefore, ( ^.llD will follow once we show 
D_D{x)R{x) = {D{x) + C) diag (e 



-ri 



X — I] 



(2.11) 
(2.12) 
(2.13) 
(2.14) 
(2.15) 

R{x)) + Vi,{x)Q. (2.16) 



(We used ( |2.4D to simplify the rhs.) 

To prove ( |2.16 ), consider its nth component. Canceling terms — exp(— r„)d„(2;)i?„ 
on the Ihs and rhs, it reads 



e'""(i„(2;)i?„(x) = dn{x)e ''"Rn{x - i) 



N 



+ Cnoer'^ {Rj{x -i)- Rj{x)) + Vb{x). 



(2.17) 



In view of ( |2.2| ) and ( p.8D , this can be rewritten as 
e'""(d„(2;)i?„(x) - dn{x - i)i?„(x - i)) 



N 



J](a,(e--K"-e-^^]) + l) {R,{. 



(2.18) 



Now from the definition (1.5) of the Cauchy matrix we see that this amounts to 



TV 



dn{x)Rn{x) - dn{x - i)Rn{x - i) = '^CnjiRjix - i) - Rjix)). 



(2.19) 



By virtue of (p.4|), this is clearly true, so (2.1(:) follows. 
Proposition 2.2. Introducing the A AO 

B = -i{Ti + Vb{x)), 
one has 

W{x,p) = {BW){x,p)+iePW{x,p). 



(2.20) 
(2.21) 
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Proof. From (|2.9| ) we obtain 

N 

Using ( |2.11| ), this can be rewritten as 



n=l 



N „ . 



(2-23) 



+ e^^PS(x - i) - gPe'^P^ fp"i^e-rl^ " Vb{x){e'''P - W{x,p)). 

n=l 



Prom the definition ( |2.8D of Vf, we now see that we can cancel three terms on the rhs. 
Then we are left with 

iW(x,p) = -eP(W(x,p) - e^^P) 

(2.24) 

+ eP {e-PW{x - i,p) - e"P) + Vb{x)W{x,p), 

which amounts to ( p.2lD . I 

Since the wave function satisfies the AAE ([1.2|), and iB equals A — Va{x)T^i (cf. ( |1.1D ), 
an alternative formula for the time derivative reads 

W{x,p) = iVa{x)W{x + i,p)- ie-PW{x,p). (2.25) 

We now proceed with the time derivatives of the potentials. 

Proposition 2.3. One has 

tix) = iil-Vaix)), (2.26) 

Vfe(x) = i{Vaix) - Vaix - i)). (2.27) 



Proof. In view of (p.8|), the time derivative (2.27) is immediate from ( p.26| ). To prove 
( 2.26| ), we first use the relation 

W(x,r„) = e^"'-"d„(x)i?„(x), (2.28) 
and the W-AAE to deduce 

e-'-Rnix + e''-Va{x)Rn{x + + {Vb{x) - e'" - e-''-) i?„(x) = 0. (2.29) 
Therefore, ( p. Ill ) can be rewritten as 

iRn{x) = e''"{Rn{x)-Va{x)Rn{x + i)), n=l,...,N. (2.30) 



Taking now the sum of these N equations and using ( |2.5D - (|2.7| ), we obtain 

itix) = A(x) - 1 - -^^[A(x + i)-l]= Vaix) - 1, (2.31) 
as asserted. I 
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In order to obtain the time derivative of Va(x), and also for later purposes, it is expedient 
to introduce one more quantity, namely, the r-function 

r(r,/x;x,t) = \1n + C{r)D{r, f,;x,t)-^\ . (2.32) 

In view of the alternative representation I(C.30) for A(x), we readily obtain 

A(x) = r(x - i)/r(x), (2.33) 

Vaix) =T{x + i)T{x-i)/T{xf. (2.34) 

We are now prepared for our last proposition. 
Proposition 2.4. One has 

f{x)/T{x) = -iT.{x), (2.35) 

Va{x) = iVa{x){Vk{x + Vk{x)). (2.36) 

Proof. Clearly, (|2j3|) follows from (|]3|), (|2j|) and (^). To prove ( p^Bj ), we first use 
( 2.32| ) and Leibniz' rule to write 



N 

T[X) = —I 

n=l 



Y^Tnix). (2.37) 



Here, t„(x) denotes the determinant of the matrix obtained from In + CD{x)^^ when 
the nth column is replaced by (e^" — e"*"") dn{x)^^{Cin, ■ ■ ■ , CnuY- In the determinant 
quotient r„(x)/T(x) we now multiply both matrices from the right by D{x). Then we 
obtain |il„(x)j/|L'(x) + C|, with r2„(x) the matrix defined in the paragraph of Part I 
containing (C.6). From I(C.7) we then deduce (2.35). I 

We now turn to some immediate consequences of the above formulas. First, we have 
from (|L1|), (|]2^) and ( p6| ) 



(2.38) 



A = Va{x)T^i + Vb{x) = iVa{x){Vb{x + i) - Vb{x))T.i 
+ i{Va{x) -V,{x-i)) = i[Va{x)T.i,T, + Vi,{x)]. 

Recalhng ( ^3oD , this can be rewritten Lax type equation 

A=[B,A]. (2.39) 
Second, from ( ^ ), ( p6D and p^SiD we deduce 

f{x)T{x) - f{xf = T{xf - t{x + i)r(x - i). (2.40) 
Third, introducing 

^'(r, fi; x,t) = i ln(A(r, fj,]x,t)) = i ln(r(r, ij,;x — i, t)/T{r, /j,; x,t)), (2-41) 
we readily obtain 

^{x) = i exp{i[^{x + i) - ^{x)]) - i exp{i[^{x) - ^{x -i)]). (2.42) 
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The nonlocal evolution equation ( 2.42| ) is the Toda type equation announced in the 
Introduction, cf. (|1.7|). Since Lemma 2.2 in Part I yields 



lim A(x) = 1, lim A(x) = exp ( 2 > ) , (2.43) 

Re X'— >oo Re x— >— oo 




it follows from ( p. 41 ) that we have 

N 

^(x) = 2i 

Re x— >oo ^ ' Re x^—oo 



lim ^{x)=0, lim ^'(x) = 2i r„, (mod 27r). (2.44) 



The multi-valuedness indicated here is inevitably present when we let x and t vary 
over C. Indeed, r(x,t) (|3|) has zeros in general (as well as poles whenever /i(x) is 
non-constant). Therefore, ^'(x,t) has logarithmic branch points. Note, however, that this 
multi-valuedness is of no consequence in ( p. 421) . 

We continue by studying reality restrictions. Specifically, we ask first: Can one choose 
the spectral data (r, /i) such that \I'(r, /i; x, t) is real-valued for real x and il 

This question can be answered in the affirmative by using the results of Appendix D 
in Part I. We showed there that the AAO A ( |l.l| ) is formally self-adjoint on L^(M, dx) 
whenever ri, . . . , r^v are purely imaginary and the functions iexp(— r„)/i„(x), n = 1, . . . , A'", 
are real- valued for real x. Along the way, we obtained as another consequence of these 
restrictions the relation 



A(x) = l/A(x), xgC, (2.45) 

cf. I(D.17). Imposing these restrictions and choosing t real from now on, the time- 
dependent factors exp(it[exp(rn) — exp(— r„)]) belong to (0, oo), so we deduce 



A(x, t) = 1/A(x, t), X, i € R. (2.46) 

Now A(x,t) is meromorphic in x and t. For x, t real, ( |2.46| ) entails that X{x,t) is a 
phase factor, so in particular no zeros or poles occur. As a consequence, we need only fix 
the logarithm branch in ( p.4lD by requiring 

lim ^'(x,t) = 0, t G M, (2.47) 

to obtain a real-valued, real-analytic function 

^-iM^^R, (x,t) ^(x,t). (2.48) 

Keeping t real, its (multi-valued) continuation to complex x satisfies 



q>{t,x) = ^{t,x), (mod27r). (2.49) 

Let us now consider the characteristics of the solutions to ( 2.42| ) with the above reality 
restrictions in force, i.e., 

Rer„ = 0, n = l,...,A^, (2.50) 
Re(e-''"/u„(x)) = 0, x G M, n = l,...,A^. (2.51) 
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Taking first = 1, and choosing ri = in'^ , £ (0, tt), we get from d^ll) and (|ll|) 

T{x,t) = 1 + (2isinK+)"^exp(-2K+x + 2fsinfi;+)/ii(x)"S (2.52) 

t) = ^ln( ^^^~'^y^^-) + exp(-2.+. + 2t sin .+) \ ^^^^^^ 

\2ie~^'^ sin K+ + e"**^ exp(— 2«;+x + 2t sin k+) / 

Since the function ie~^'^^ fj,i{x) is real- valued for real x, this indeed yields a real- valued 
function for S M^. But whenever fj.i{x) is non-constant, ^{x,t) cannot be viewed as 
a 1-soliton solution. Indeed, in that case it is not of the traveling wave form f{x — vt). 
Choosing however 

2iie"*'""^/ii(x)sinK+ = exp(-2K+a+), a+ G R, (2.54) 

we do get a function that is not only of the form f{x — vt), but also of the kink type, in 
the sense that its x-derivative is positive and exponentially localized around its maximum 
at X = a+ + vt: 

, , , 2k+ sin , , sin k 

^{x,t) = r-rvT, v{k) = . 2.55 

^ ' ' cosK+ + cosh2K+(x-a+ -7;(K+)t) ' k ^ ' 

It is readily checked that the velocity function w(k) decreases monotonically from 1 to 
as K goes from to vr. Observe also that one has 

lim ^'(x,t) = -2k+, teM. (2.56) 



Consider next the choice r\= ik — in, k G (0, vt). Then we obtain 

t{x, t) = 1 — {2i sin k~)~^ exp(— 2k~x — 2t sin K^)fii{x)^^ , (2.57) 

^, , , / -2ie^*'^ ui(x) sin K~ + e^'^ exp(— 2k"x — 2t sin k~) \ 

= ^ 111 ^. . _^ V ^ _ TT • (2-58) 

\—2te ^'^ fj,i[x)smK + e ^'^ exp(— 2k x — 2rsmK j / 

Once again, we need fJ-i{x) to be constant for this function to have 1-soliton character- 
istics. Setting 

- 2ie"*''~^i(x) siuK" = exp(-2K"a"), a" € M, (2.59) 

we now obtain 

, 2K~sinK~ 
x,t = 7 r— VT- 2.60 

cosK +cosh2K [X — a + v[k )t) 

Thus, ^(x,t) is a 1-kink solution moving to the left with speed in the interval (0,1). 
Clearly, it satisfies 

lim ^'(x,i) = -2k", teR. (2.61) 

X—^ — 00 

Proceeding with the arbitrary- case, we can clearly make choices for Unix) correspond- 
ing to (|2.54| ) and (|2.59| ), depending on whether r„ belongs to i(0,7r) or «(— 7r,0). Doing 
so, we obtain real-valued, real-analytic solutions. In Section 6 we study these solutions in 
detail, demonstrating in particular that they may be viewed as A^-soliton solutions. 
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3 Reflectionless self-adjoint Schrodinger operators 

In this section we clarify the relation of the above AAOs and their reflectionless eigen- 
functions to the reflectionless self-adjoint Schrodinger operators considered in the 1ST 
framework |11, 12, |l^. Thus, we start from a continuous real-valued potential Vh{x) with 
decay at ±oo given by 

Vh{x) = O {\x\-^) , x^±oo, (3.1) 

with a suitable positive exponent e. Then it is clear that the Schrodinger operator on 
-L^(IR, dx) given by 

{Hf)ix)^-f"{x) + VHix)f{x), (3.2) 

is self-adjoint on the natural domain of the free Hamiltonian —cP/dx'^. 

The reflectionless operators can now be characterized by the existence of an H-eigen- 
function 

{HWh){x,p)=p^Wh{x,p), (3.3) 
with asymptotics 

I exp(ixp), X ^ oo, , , 

Wh(x,p)~ 7/ , (3-4) 
I ohCpJ exp(ixpj, X ^ — oo. 

The 1ST framework yields a complete classification of such operators: The function anip) 
is of the form 

N 

aHVP) = — —. — , < < • • • < Ki, (3.5) 

P + lKn 

n=l 

and for each such S'-matrix there exists an A^-dimensional family of potentials parametrized 

by normalization coefficients vi,...,!/^ G (0, oo). The eigenfunction Wh{x^p) yields 
bound states 

</>«(•) = W//(-,^^n), n=l,...,A^, (3.6) 
satisfying 

|</.„(x)|2(ix = IK, n = l,...,A^. (3.7) 



Thus H has continuous spectrum [0, oo) and discrete spectrum —k\,..., — k^. The bound 
states are pairwise orthogonal, and the improper eigenfunction Wh{x-,p),p € M, gives rise 
to the (Schwartz) kernel of an isometry from L-^(]R, dp) onto the orthocomplement of the 
bound state subspace. 

As we will show in Part III, the state of affairs concerning Hilbert space properties of 
the above AAOs and their reflectionless eigenfunctions deviates from this "Schrodinger 
scenario" in several ways. In this section, however, we only aim to show how the self- 
adjoint reflectionless operators H and their eigenfunctions W//(x,p) arise as limits of our 
AAOs and their reflectionless eigenfunctions. 
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Turning to the details, we fix 

{k,u) e {0,oof^ , < Kjv < ■•• < Ki. (3.8) 
Then we choose 

rn = ij3Kn, /3 G (0,7r/Ki), n = l,...,N, (3.9) 
and constant multiphers 

^Jin{x) = -ie'" / dvn, n = l,...,iV. (3.10) 



This entails that the A AO A ( |1 . 1[ ) is formally self-adjoint, cf. Theorem D.l in Part I. 
Moreover, since all numbers ri, . . . , r^v have imaginary part in (0, vr), we may invoke Theo- 
rem 3.3 in Part I. It entails in particular that A can be rewritten as 

A = Sl-2, (3.11) 

where 

5+ = T,/2 + F(x)T_,/2, (3.12) 

N 

V{x) = Y,iRn{x - i/2) - Rn{x)) + 1. (3.13) 

n=l 

Furthermore, the wave function is an S'^-eigenfunction: 

{S+W)ix,p) = [fPl^ + >V(x,p). (3.14) 

We now introduce a scaling x P^^x, p ^ Pp in the various quantities at hand, and 
study the /? — > limit. Specifically, introducing first 

Rp^^ix) = p-^Rn {p-^x) , (3.15) 

the system (2.4) with x P~^x can be rewritten as 

exp(i/3«;„ -|- 2k„x) 



N 



P 

~^ ^ exp(z/5K„) — ex.p{—ipKj)^^'-'^^^ ^ 1, . . . , iV. 



(3.16) 



From this one readily deduces that Rp^n{x) is holomorphic at /? = 0, with a limit 
jm.^Rp^n{x) = R^{x), n = l,...,iV, (3.17) 

that solves the system 

N 



exp(2K„a;) jj 



R^ix) 



+ . \. Rf{x) = 1, n = l,...,N. (3.18) 
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Secondly, setting 

Wp{x,p)^W{p-'x,l3p), (3.19) 
it follows that Wf3{x,p) is holomorphic at /3 = as well, with a limit 

hm Wf3{x,p) = 1 - E ^ ^h{x,p). (3.20) 



\ .7 = 1 



'3 



Thirdly, introducing 

5^,+ = exp(-i/?9,/2) + Vp{x) exp(i/3a,/2), (3.21) 

with 



AT 



Vp{x) = /3 - W2) - i?/3,n(x)) + 1, (3.22) 



n=l 



we clearly have 

S^,+W^{x,p) = + W^(x,p). (3.23) 

Now from ( |3l7| ) we deduce 

F;3(x) = 1 + {p/2fVH{x) + O {/3^) , /? ^ 0, (3.24) 

where 

N 

VH{x) = -2iY,d.R^{x). (3.25) 

n=l 

As a consequence, we obtain 

5/3,+ = 2 + (/3/2)2f + O (/?3) , /3 ^ 0, (3.26) 

with 

H = -dl + Vnix), (3.27) 

and 

(//>Vh)(x,p) =/>v^(x,^5). (3.28) 

Proceeding in this way, we have actually obtained all of the reflectionless self-adjoint 

Schrodinger operators delineated at the beginning of this section. To transform them to 

a more familiar form, we need only invoke ( |3.20| ) and the system ( ^.18 ) to write 



-B^{x), n = l,...,iV. 
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Thus Vh{x) ( ^.251) can be rewritten as 



N 



Vh{x) = 2 ^ und. {e-^-^cl)n{x)) . (3.30) 

n=l 

This is the weh-known formula expressing the reflectionless potentials in terms of their 
bound states (/>i(x), . . . , (j)N{x). 

To conclude this section, we add three remarks. First, we point out that the asymptotics 

I exp(ixp), a; ^ oo, , ^ 

W{x,p)^{ p . (3.31) 

N 

<P) = n -V (3.32) 

n=l 

derived in Theorem 2.3 of Part I yields the function anip) ( |3.5[) when one substitutes 
p (3p, rn i(3nn, and takes /? — > 0. Second, we observe that we could also have used 
the AAO 

Ap ^ - 2 (3.33) 

to arrive at the Schrodinger operator H, as will be clear from the above. 

Finally, just as the operator H may be viewed as a reduced Hamiltonian for a Galilei- 
invariant two-particle system, one may view the AAO 5/3^+ (or A^) as the reduced Hamil- 
tonian for a Poincare-invariant two-particle system. In this scenario the limit /3 — > 
amounts to the nonrelativistic limit, cf. Ref. |14]. 



4 Reflectionless self-adjoint Jacobi operators 
and Toda lattice solitons 

In this section we show how the reflectionless self-adjoint Jacobi operators and Toda lattice 
solitons arise via analytic continuation and x-discretization. For comparison purposes, 
especially useful are the monographs Refs. |l^] and Flaschka's paper Ref. |17]. For the 
Jacobi operator 

(J/)(n) = a(n - l)/(n - 1) + a(n)/(n + 1) + 6(n)/(n) (4.1) 

on i^(Z) one requires decay 

a(n) = 1/2 + (|n|-^°) , b{n) = O {Inl'^^) , n ^ ±oo, (4.2) 

for suitable positive exponents Bq, ei,. Requiring in addition that a{n), h[n) be real, it is 
clear that J is a bounded self-adjoint operator on /^(Z). 

With these requirements in effect, the reflectionless Jacobi operators are characterized 
by the existence of a J-eigenfunction Wj(n,p) satisfying 



{jy^j){n,p) = cos{p)Wj(n,p), 



(4.3) 
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with asymptotics 

Wj{n,p) ~ exp(mp), n — >• oo, (4-4) 

Wj{n,p) ~ aj(p) exp(mp), n — > —oo. (4.5) 
From the 1ST formaUsm it follows that the function aj{p) is of the form 

yr Sin[(p - iKp/2] ^ COs[{p-iK-)/2] 

where < /t^^ < • • • < /tf , (5 = +, — . For each aj{p) there exists an {Nj^ + A^_)- 



' ^Af^ ' = +5 — • The operator J has + bound states, given by 



(4.7) 



dimensional family of Jacobi operators, parametrized by positive normalization coefficients 
1^1) • • • ) ^Ns^ 

(t>+{n) = Wj{n,iK+), j = l,...,N+, 
(j3Y{n) = Wj{n,TT + iKl), / = l,...,iV_ 

and normalized as 

^|0i(n)|' = l/i.f, k = l,...,Ns, <5 = +,-. (4.8) 

ngZ 

The bound state energies equal 5 cosh(Kj), j = 1, . . . , N^, 5 = +,—. 
Fixing the above spectral data, we choose 

iV = iV+ + iV_, (4.9) 

and let 

rj=e'^K+, i = l,...,A^+, rjv++/ = e^%- - z^, Z = l,...,iV„. (4.10) 

Here, we choose r] G (Ojr/^), with r]s € (0, 7r/2) satisfying K^sinr/s < tt, 5 = +,— . Then 
the requirements ( |l.3D and (|l.4| ) are clearly met. We also choose constant multipliers 

^j(x) = exp(rj)/z^+, j = 1, . . . , iV+, 
^Af++/(2;) = exp(rAr^+i)/z^;", / = 1, . . . , iV_. 

Next, we substitute 

X ^ ie'^'^n, p^-ie^'^p, n G Z, r] e {0,7]s), (4.12) 

in the above quantities, and study the t/ ^ limit. To this end, it is convenient (both for 
notation and for comparison purposes) to introduce further parameters 

^^fexp(-4), ,.^/"^ J = l,---,N^, 

' \ -exp(-/tT_^^), ' \ "i-N^^ j = N+ + l,...,N. 

Employing these parameters, one easily checks 

lim D (ie-^^n) = diag (1/^2"+ Vi, . . . , 1/4"+ V^v) = L'^(n), (4.14) 
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lim C,k = z = C-L/zk, j,k = l,...,N. (4.15) 

rj^O 1 - ZjZk 



Using Zj G (—1, 1), the matrix C is readily seen to be positive. (This can be deduced from 
Cauchy's identity. Alternatively, following Flaschka |17], one need only write 1/(1 — ZjZk) 
as a geometric series to verify positivity.) 
As a consequence, the matrix 

L\n) = \\raCD{ie-''^ny^ , n e Z, (4.16) 
is given by 

L-^(n),fc = C/fczf z.fc, j,A; = l,...,iV, n G Z. (4.17) 

Therefore, L'^{n) is similar to a positive matrix, so Iat + L'^{n) is invertible. From this 
and the limit ( |4.14[ ) it readily follows that 

lim Rj {ie-'^^n) =Rj{n) (4.18) 



exists for all n E Z. Since Rj (n) solves a system 



TV 



1 •''J 



+ Y,—^Ri{n) = l, j = l,...,iV, (4.19) 



k=l - ^^^'^ 



with real coefficients, it is real, too. 
Next, one easily verifies 

liiaD {ie-'''n±i) = D-^{n±i), n G Z. (4.20) 

From this it is not hard to deduce 

lirnRj {ie-''>n±i) = Rj{n±i), j = l,...,N, n E Z. (4.21) 

(Indeed, the limit of the system (^3) in small complex neighborhoods of x = in ± i can 
be controlled by exploiting ( |4.19D with n — > n ± 1.) 

Using the above, the remaining pertinent limits can be easily obtained. First, recalling 
(|23^)-(|23^) (with t = 0), one gets 

lim T (ie-''^n) = \1n + L-^{n) I = T-^(n), (4.22) 

rj— >0 

lim A (ie-'^'n) = r^n - l)/r-^(n) = \\n), (4.23) 

lim Va (ie~*^n) = X\n)/X\n + I) = {2a{n)f . (4.24) 

r;— »0 

This defines a(n) E (0,oo), since d^) entails T-^(n) > 0. Secondly, (U) yields 

N 

lim Ffc {ie-'^^n) = ^ (i?/ (n - 1) - i?/ (n)) = 2b{n). (4.25) 
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Since Rj{n) is real, b{n) is real, too. 
Thirdly, consider the wave function 

W^{x,p) = W {ie~''^x, -ie'ip) . (4.26) 

Its limit 

/ ^ R-^(n) \ 

lim W,(n,p) = e^"P ^ ~ ^ e-i - z^ I ^ ^''^^'P^' ^ ^' ^^'^'^^ 

is once again immediate from the above. Now when we fix n G Z, the AAE ( |1.2D yields 

Wrj (n - e'^,p) + K (ie-^^n) {n + e^^p) 
+ {Vb {ie-^'^n) - 2 cos (e*'?p)) Wr,{n,p) = 0. 

Taking to 0, we deduce that >V"^(n,p) satisfies the discrete difference equation 

yV-^{n -l,p) + iainfW^n + l,p) + 2b{n)W\n,p) = 2cos{p)W\n,p). (4.29) 

This is not yet of the Jacobi form (O), cf. (^J]), and we now proceed to explain how the 
connection is to be made. 

Our reasoning involves the |n| — > oo asymptotics of the above quantities. This asymp- 
totics does not follow from our previous work. For one thing, we have taken a limit that 
goes beyond the parameter regime studied in Part I. But even when one ignores this, the 
pertinent asymptotics is quite different from the one in Part I. Indeed, here we need the 
asymptotics in the direction of the shifts n ^ n it 1, whereas in the A AO context we are 
dealing with the asymptotics Re x — > ±00, which is orthogonal to the shifts x ^ x ± i 
in (PI) . 

It is however a simple matter to obtain the desired n ^ 00 asymptotics directly. Indeed, 
from (|1|) and ( p^) one readily deduces 

Rj{n) = 0{ex.p{-pn)), n ^ 00, j = l,...,iV, /j = 2min(K+_^, k^_). (4.30) 



Likewise, from (4.17) one has 

L-^{n)jk = 0{exp{-pn)), n ^ 00, j,k = l,...,N. (4.31) 
Then (|]22|) yields 

n —> 00, (4.32) 
so from (|423|) and (g^ ) one obtains 

A-^(n) = l + 0(exp(-pn)), n ^ 00, (4.33) 

2a(n) = 1 + 0(exp(-pn)), n ^ 00. (4.34) 

Hence the infinite product 

00 

n(n) = JJ (2a(m))-^ (4.35) 
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converges, and one has 
n(n) = A-^(n)-i/2_ 

When we now renormahze yV'^{n,p) by introducing 

then it fohows from ( [4. 291) that we have 

a{n - l)W^{n - l,p) + a{n)W^ (n + l,p) + b{n)W-^{n,p) = cos{p)W-^ {n,p), (4.38) 



(4.36) 



(4.37) 



which is of the Jacobi form ([4.3|). It is therefore clear from the above that the n 
asymptotics is in accord with (|]2|) and (Q: From (|430| ) and (|]2|) one gets 



oo 



b{n) = 0(exp(— /on)), 



n oo, 



(4.39) 



which, together with ( 4.34 ), agrees with ( [4.2| ) for n — > oo. Moreover, (|4.30D , ( ^2^ and 
(|434|)-(|43^) entail (gj). 

It remains to show that the asymptotics for n —oo works out as announced. This 
involves a little more work. Let us first note that (4.19) entails that {n) has a finite 
limit i?"^(— oo) for n — > — oo, satisfying 



N 



fc=i ^ ^^^^ 



i = l,...,iV. 



Using Cramer's rule, this can be improved to 



rJ (n) = R^- (-oo) + 0(exp(pn)), 



n 



-oo. 



l,...,iV. 



Therefore, ( [4. 25] ) yields 

h{n) = 0(exp(/?n)), n — > — oo. 
Next, from ( ^p2^ , ( 14221) and (p^ ) we have 



(4.40) 



(4.41) 



(4.42) 



A^(n) 



diag (z^ 



2n 



-2n\ 



diag [zl 



■2n 



-2n\ 



(4.43) 



N 



llz~^ + 0{exp{pn)), 



n 



-oo. 



k=l 



Therefore, we obtain from (4.24) 
l/2 + 0(exp(/>n)). 



an 



n 



-oo. 



(4.44) 



Summarizing, the asymptotics for \n\ 
with (O). 



oo of the coefficients a{n) and 6(n) agrees 
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We are left with determining the n ^ —oo asymptotics of the wave function. A quick 
way to obtain this is to invoke the alternative representation for yV(x,p) from Theorem C.3 
of Part I. It entails that yV'^{n,p) ( ^2% ) can also be written 

W^in p) = ^^n, \D'in) + C^^^-E{z^\...,z-,')A^ip)\ 

^ '"^^ |l)» + C^diag(zr\...,%^)| ' ^ ^ ^ 

with 

A^ip) = diag {S-^{z,;p), d'^{zN;p)) , (4.46) 
SJ^z; p) = l- ~ ^ = '^~'''-^~\ (4.47) 
Thus we have 

TV 

Vy^(n,p) ~ e*"P n^-oo. (4.48) 

From ( [4.37| ) and ( [4. 431) we now obtain 

^ e-*P - z"^ 

W^(n,p) ^e'"P]T|z,-| ' , n^-oo. (4.49) 

XX e ^P — Zj 

j=i ■> 

Substituting ( [4.13| ), this yields (pl5|)-(p^, as announced. 

With the above r] ^ limits at our disposal, it is straightforward to calculate the limits 
of the time-dependent quantities and time derivatives in Section 2. There is however one 
crucial change to be made before doing so: We should replace t by it so as to obtain the 
pertinent real-valued Toda lattice quantities for r/ — > 0. 

Indeed, doing so in the formula ( 2.10[) that defines the time-dependence, one obtains 

limL'(ie-^''n,«t) = diag(l/z2"+iz/i(t),...,l/4"+ii/^(t)) =D-^{n,t), (4.50) 

i^j (t) = exp [-t (zj - ^7^) ) , j = 1, . . . , A^. (4.51) 

The time-dependent version of ( [4.16 ) then yields the limit 

L-^(n,t)jfc = C/fczfi.fc(t), j,fe = l,...,Af, neZ, (4.52) 
and ( 4. 22] ) is generalized to 

lim T (ie~*'^n, it) = IIat + L^(n, t)| =T\n,t). (4.53) 

Then ( p. 40 ) entails Hirota's formula [|l8| 

Likewise, (pIS^ ) and ( ^^271) , together with ( |424| ) and (g3|), yield 

a{n,t) = a{n,t)[b{n,t) - h{n + l,t)], (4.55) 
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h{n,t) = 2a{n - l,tf - 2a{n,tf 



(4.56) 



Comparing with Flaschka's paper Ref. [17|, we see that ( 4.55| ), (|4.56 ) coincide with his 
Eq. (2.3), up to signs. This sign difference arises from our different sign convention for 
the time-dependence in ( [4.51| ). (Our convention agrees with Refs. ||T^, 16 1.) Denoting the 
quantities he uses in his Section 3 with a superscript F, one gets (again up to irrelevant 
conventions) 



4)' 



j = l,...,N, 



-1 AFn 



n £ Z, 



(4.57) 
(4.58) 
(4.59) 



(Compare his system Eq. (3.3) to ( 4.19| ) and note his Eq. (3.7) to check these correspon- 
dences.) 

Various formulas in Refs. |15, |l^] can also be obtained as — > limits. In this con- 
nection we mention in particular the Toda/Kac-van Moerbeke account in Section 3.8 of 
Toda's monograph Ref. [^]: The pertinent formulas can be readily derived by taking the 
r] ^ limit of results that can be found in Section 3 of Part I. 

Furthermore, we point out that van Diejen recently obtained the analog of the formu- 
la ( 4.45| ) for the Jacobi wave function \V'^{n,p), cf. Ref. ||l^. His paper also contains further 
results of interest, and information on recent literature dealing with Jacobi operators and 
the Toda/Kac-van Moerbeke correspondence. 



5 Parametrization via relativistic Calogero— Moser systems 

Returning to the general setting of Section 2, we proceed to detail the connection with 
the IIrei(''" = "^f^) systems from Ref. |^]. The connection to these finite-dimensional 
soliton systems hinges on a suitable reparametrization of the Cauchy matrix ( |1.5D and 
the multipliers in the diagonal matrix D{r,fj,;x) (p.l|). To ease the notation we assume 
from now on that the numbers ri, . . . , r^v are ordered such that ^i, . . . , r7v+ have imaginary 
parts in (0, vr) and rAr_Ar_+i, . . . , r^r in (— vr, 0), with 

iV+ G {0,1,... ,Af}, N. = N-N+. (5.1) 

(Since the wave function W(r, x,p) and AAO A{r, /x) are invariant under permutations 
on the data (r, //), this does not give rise to a loss of generality.) 

It now turns out that the numbers ri, . . . ,riy^ and r7v-7V_+i, ■ ■ ■ i^N can be traded 
for the (complex) positions qf , . . . and , . . . , q^_ of the particles and antiparticles. 



resp., in the Ilrei system, in such a way that the reparametrized Cauchy matrix C(r) (1.5) 
and the Cauchy matrix C in the Ilj-ei Lax matrix are closely related. This relation then 
suggests a reparametrization of the multipliers in D{x) such that the matrix CD{x)~^ 
may be reinterpreted as the IIrei(T = '^/'^) Lax matrix evaluated in x-dependent points of 
the (complexified) IIrei(''" = t^/"^) phase space. 

More precisely, the latter identification holds true up to diagonal similarity transforma- 
tions. Since we are dealing with determinants and spectra, such "gauge transformations" 
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are immaterial. In particular, the definition Eq. (2.70) of the Lax matrix in Ref. 
amounts to a gauge choice that facilitates its spectral analysis, but in the present context 
another gauge is more convenient. 

Specifically, here we work with the Lax matrix 

£((?, 6) ^ C{q+,q-)V{q+, q-,9+, 9~). (5.2) 

The Cauchy matrix C is defined by 

C,-fc = l/cosh[(g+-g+)/2], (5.3) 

CN++i,N++m = l/cosh[{qf - g„)/2], (5.4) 

CN++l,k = -i/ sinh[(g~ - q'l)/2], (5.5) 

Cj,N++m = i/ sinh[(g+ - q~)/2]. (5.6) 

Here and from now on, the indices j, k take values 1, . . . , A^+, whereas the indices I, m 
take values 1, . . . , The diagonal matrix T) is defined by 



V ^ diag (^exp((?+)y+, . . . , exp(e+jy+^, exp(^r)Ff , . . . , exp(e^_)y^_ j . (5.7) 

The quantities 9^ are the generalized momenta corresponding to the positions q^, and the 
"potentials" are given by 

H |coth[(g;-g+)/2]| H |tanh[(g;-g-)/2]|, (5.8) 
n |coth[(g--gJ/2]| II |tanh[(g--g+)/2]|. (5.9) 

l<m<N-,rri=/=l l<j<N+ 

(The moduli we are choosing here preclude analyticity in q, but they enable us to steer 
clear of multi-valuedness issues. Such issues are important in other contexts, but here 
they would give rise to unnecessary complications.) 

Comparing (^.2|)~ (|5.9| ) with real q's and ^'s to the Lax matrix L given by Eq. (2.70) in 
Ref. Q, we see that C amounts to a diagonal similarity transform, as announced. More 
precisely, we should substitute r = PfJ.g/2 = 7r/2, ^ = 1, /3 = 1, and xf^ q^, pf^ ^ 9f^ 
in loc. cit., and choose distinct qf, . . . cf. also Eqs. (1.3) and (1.4) in Ref. [^. (To 

avoid confusion, we should add that the Cauchy matrix C (|5.3D-(^^) is slightly different 
from the matrix we refer to as Cauchy matrix in Ref. The latter equals £C6, with 



£ = diag (^exp(-g+/2), . . . , exp(-g+^/2), exp(-gi /2), . . . , exp(-g^_/2) j , (5.10) 

cf. also Eq. (Bl) in Ref. ||.) 

After these preliminaries, we turn to the reparametrizations of C(r) and the multipliers 
in D{r,fi;x). To this end we first introduce parameters 

= -irj, (5-11) 
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a, = -irN++i + TT, 
which are convenient in their own right. Clearly, one has 



Re € (0, vr) 



n = 1, 



.Ns, 



+, 



and the Cauchy matrix (1^) can be rewritten as 
Cjk = 



1 1 

- exp(-ia+/2 + ia^/2)—-——-^—-, 

2 ^ sm[(aj + a^)/2] 



i _ _ 1 

CN++i,N++m = ^ exp(-ia; /2 + ia^/2)—-— — — -— , 

CN++i,k = -I exp(-ia,-/2 + ia+/2) 

2 cos[{ai +al)/2] 



1 



1 



Cj,N++m. = exp(-ia+/2 + ia^/2) _ 

2 cos[(a7 + am)/2J 



We also note that the restrictions (1.4) amount to 



a| + a, 7^ vr. 



(5.12) 

(5.13) 

(5.14) 
(5.15) 
(5.16) 
(5.17) 

(5.18) 



(Recah j,k G {1, . . . , A^+} and /,m € {1, . . • ,iV-}.) 

Next, we relate the quantities and to the above positions and g^", resp. To 
this end, we begin by pointing out that the map a ^ z = cot(a/2) yields an injection of 
the strip Re a G (0, vr) onto the half plane Re z > 0. (Even though this assertion is not 
immediate, it is straightforward to check.) Therefore, the map 



F : {Re a e (0,7r)} ^ {Im g G (-^/2, 7r/2)}, 



a 



ln(cot(a/2)), (5.19) 



with F{tt/2) = 0, is holomorphic and has a holomorphic inverse. It is easily seen that the 
resulting relation 



e"9 = tan(a/2), q = 4^ a = Tr/2, 
implies 

cosh g = 1/ sin a, sinh(7 = cota, tanhg = cosa. 
We now set 

qi = 5F{ai), n = l,...,Ns, 5 = +,-. 
Thus we have 

Img^ G (-7r/2,^/2), n=l,...,Ns, S = +,-. 
Moreover, the restrictions ( ^.18 ) translate into 



(5.20) 
(5.21) 
(5.22) 
(5.23) 
(5.24) 
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In words, the positions qf , . . . , Qp^^tQi , ■ ■ ■ iQi\f must be distinct. 

When we combine the changes of variables r {a~^,a~) iq^,Q~), we get from the 
above 



exp(rj) = {i + sinhqj')/ coshg^, 
exp{rN^+l) = {-i + sinhg'^')/ coshg^" 



(5.25) 
(5.26) 



But when we would substitute this directly in the Cauchy matrix ( |1.5D , we would obtain 
a matrix whose relation to the Cauchy matri x C (^.3| )-( |5.6| ) is invisible. Instead, we start 
from ( |5.14| ) and use the relations ( ^.20 ) and ( 5.211) to calculate 

sin(a+ + Q+)/2 = cos(a+ /2) cos(a^ /2) tan(a+ /2) + tan(a+ /2) 



1 / sin(a+) sin(a^) 

2 \^tan(a+/2)tan(a^/2) 

cosh[(g+ - g+)/2] 



3 
1/2 



tan(a+/2) +tan(a^/2) 



(5.27) 



[cosh(g'+) cosh(g^)]V2 

(Note that the radicands have no zeros or poles in the pertinent regions, cf. ( |5.19D . The 
branch choice is then clear: We need the positive square root for real.) 

More generally, in terms of q^, q~ the Cauchy matrix C ( 5.14| )-( [5.17| ) becomes 



C = SCS~^D. 



where 



Dc = - diag y — cosh q^, . . . , — cosh q]^_^ , cosh q^ , . . . , cosh q^ 
S = diag ( exp(— ia^/2)[cosh qf]^^"^, . . . , — exp(— ia^_ /2)[cosh qj^_]^^'^ 



(5.28) 

(5.29) 
(5.30) 



We stick to the parameters , in the exponents, since this will be convenient shortly. 

We now turn to the reparametrization of the multipliers in the diagonal matrix D{x). 
In terms of a"*", a~, this matrix reads 

(5.31) 



(5.32) 



D{x) = diag yjii{x) exp(2a||^x), . . . ,^n{x) exp(2a^_x) 
cf. (|2lD-(0) and (|5]l|)-(|5]l|). Consider now the matrix 

L{x + i/2) = S-^CD{x)-^S = CDcD{x)-\ 



where we used ( 5.2^ ). Comparing it to the Lax matrix C (|5.2[) , we see that when we 
rewrite the multipliers as 



Hj{x) = -- exp(ia^) cosh(g+) V'^+7r+(x + i/2) 



-1 



= -exp{iaY)cosh{qf) [Vf7r,^{x + i/2)] 



-1 



(5.33) 
(5.34) 
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then we have 

L{x) = C{q,e)M{x), (5.35) 

with 

A4{x) = diag ^exp(— 0^ — 2a^x)7r^(x), . . . ,exp{—6]^_ — 2a]^_x)7rJ^_{x)^ . (5.36) 



The above formulas ( |5.28 )-( 5.36D encode the announced relation between the two key 



matrices C and D{x) and the IIrci(''" = '^/'^) Lax matrix from Ref. |||]. Note that the 
parametrizations ( ^.331 ) and ( |5.34| ) give rise to well-defined i-periodic meromorphic multi- 
pliers Tr^{x), with finite limits 

lim 7r^(x) = 7r^, n=l,...,Ns, 5 = +,-. (5.37) 

|Re a;|— >oo 

Admittedly, at this stage it is not clear that the relation just established is useful. At first 
sight, it merely seems a bizarre coincidence, and indeed the use of the position variables 
would have been quite inconvenient in Part I and in Section 2 of the present paper. 

As it turns out, however, the relation can be used to great advantage, not only for 
studying A^-soliton solutions (which we do in Section 6), but also for studying the above 
A AOs from the viewpoint of quantum mechanics (which we do in Part III). For both of 
these applications, we restrict attention to purely imaginary ri, . . . ,ri^, and to constant 
and positive multipliers 7^f^. The associated positions are then real, and we may and 
will choose real such that 

J _ „„„ (oS 



< = exp(^0^j, n = l,...,iV5, S = +,-. (5.38) 

(Notice that these restrictions amount to the ones at the end of Section 2.) 

With the choices just detailed in effect from now on, we obtain from each point in the 
phase space 

J7 = {((?+,..., Q- ,et,...,e], )GM2Ar|^+ . 

^ " " (5.39) 

9iv_ < ••• < 9i ,qt ^Qi ^ j = I = 1,...,A^_| 

of the IIrei(''" = T^f^) system [|| an AAO A that is formally self-adjoint, and a real-valued, 
real-analytic solution to (p..7|). We now turn to a study of the latter solutions. 

6 A close-up of the AT-soliton solutions 



The r-function ( |2.32| ) associated to a point {q, 9) in the phase space (pD can be 
rewritten as 

r(x,t) = |1^ + L(x + i/2,t)| = \1n + L{x,t)U\, (6.1) 

with 

L{x,t) = ciq-^^q'^O^ -2al[x-v{al)t],...,e^_-2a]^_[x + v{a]^_)t]\ , 

^ ' (6-2) 

v(a) = sin a. 
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U = diag ( exp(— zaj"), . . . , exp(— za^_) j . (6-3) 



(Recall (p3D-(^38D and (|2lo| ) to see this.) The corresponding solution (|2^ ) to { ^A^ 
then reads 

"^ix, t) = iln[ — — . (6.4 

V \lN + L{x,t)U\ ) 



In particular, for A'^-i- = 1, = 0, one obtains from (|5.2| )-( p^ the right-moving soliton 
/+((7+, X — v{p&^i)^ where 

rfeMj^nnf ^ -^°";'V (6.5) 

' ^ \^ 1 + e-*"+ exp(0 - 2a+x) y ^ ' 

= 2 Arctan (exp(— (7)), (6.6) 
and for = 0, = 1, the left-moving soliton f^{q^, 6~;x + v{a^)t), where 

f iq,0;x)=iln — ^T-rh (6-7) 

\ 1 + e exp(6' — 2a x) I 

a~{q) = 2 Arctan (exp(g)). (6.8) 

(We recall that we have already studied these solutions at the end of Section 2.) We proceed 
by showing that the general A^-soliton solution ^'(x,t) ( |6.4| ) has a long-time asymptotics 
that is a linear superposition of N 1-soliton solutions. 
To this end we define 

N+ 

¥'\x,t) = 5^/+ [q+,ep,x + 6A^iq+,q-)/4aj-viapt 



AT- 

1=1 

where 5 = +,— and 

A;((Z+, ^ I E - E I 1^ (coth2[(g+ - g+)/2]) , t) 
\k<j k>j) 

iV_ 



(6.9) 



(6.10) 



1=1 



(tanh2[(g+-g-)/2] 



\m<l vrOlJ 

N+ 

f^ln (tanh2[(g--g+)/2] 



(6.11) 
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Moreover, we introduce the remainder functions 
R^^\x,t) = ^{x,t) -^^^^{x,t), 6 = +, 
Setting 

we conjecture that one has 

= 0(exp(-(5tr)), 6t ^ oo, 6 
with 

(2aj\viap - v{a+)\,2a+{v{ap+vial)), 



+,■ 



r = mm 



(6.12) 



(6.13) 



(6.14) 



(6.15) 



In Section 7 of Ref. |p we proved the analog of this conjecture for the particle- Hke solutions 
to the sine-Gordon, modified KdV and KdV equations. Unfortunately, the strategy we 
followed in the latter cases does not apply here. 

On the other hand, the structure of the solution ^'(x, t) ( |6!^ ) makes it possible to supply 
a quite simple and direct proof of a weaker convergence result, expressed in the following 
proposition. 



Proposition 6.1. Fixing xq, sq € M, one has 

lim exp (-iR'^^\xo + sot,t)) =1, S = +, 



(6.16) 



In order to appreciate this proposition and to prepare for its proof, we begin by pointing 
out that the function 



Ftixo, So) = {q, 9; xq + sot - v{a'^)t) 
has the following discontinuous limiting behavior: 

0, So > via-^), 



(6.17) 



lim Ft{xo,so) 



t— >oo 



— 2a+, So < v{a'^), 

, /+(g,6';xo), So = u(a+). 



(6.18) 



(Indeed, this is clear from (|6^).) From this the t — > ±oo limits of ^'(^^(xo + sot,t) 
(cf. (|6.9|) ) are quite easily calculated. Denoting the limits by F^^\xo, sq), Proposition 6.1 
is equivalent to 



lim ^'(xo + sot,t) = -F(^^(xo,so) + 27rA;, 6 = +,-, k £ Z. 

St^oo 



(6.19) 



A priori, the integer k depends on 6, xq and sq, though it is undoubtedly true that one 
has k = (recall the branch choice ( p.47| )). But we cannot rigorously deduce this, since 
( 6.16| ) is a pointwise limit. 
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Next, we observe that the i-dependence of ^'(xo+sot, t) is carried by factors exp(tAi(so)) 
in the Lax matrix, with 

X(s)^l -^"tis-v{at)), i = l,...,iV+, 

' \ -2a-_^^{s + v{a-_^^)), z = iV+ + 1, . . . , iV, 

cf. (|6.2|). Furthermore, all principal minors of the matrices LU^^ occurring in ( |6.4| ) can 
be readily calculated via Cauchy's identity, yielding non-zero numbers. 

We now turn to a lemma in which the state of affairs just sketched is studied in a 
somewhat more general setting. For M an x matrix, denote by Mj the j x j matrix 
obtained from M by deleting the rows and columns j + 1, . . . , N . Now we define a set A4 
of matrices by 

M = {M £ Mn{C) I \Mj\ / 0, j = l,...,Af}. (6.21) 

Lemma 6.2. Let M+,M- G M and let 

Q{t) = \1n + M+e'^\/ \1n + M-e'^\, (6.22) 

where 

D = dmg{di,...,dN), di, . . . , c^at G M. (6.23) 
Assuming 

G (0,oo), dn+i,...,dN e {-00,0), n€{0,...,N}, (6.24) 

one has 

limQ(f) = |M+|/|M-|, (6.25) 

I— >oo 

with \Mq\ = 1. Assuming next 



di, . . . G (0, oo), dn+i = 0, 

dn+2, ...,dN G (-00, 0), n G {0, . . . , A^ - 1}, 



(6.26) 



one has 

hm Qit) = (|M+| + |M„VI)I/(|M-| + IM-+J). (6.27) 
Proof of Lemma 6.2. We may write 

Q{t)=r]+{t)/7]-{t), (6.28) 

where 

r]\t) = |diag {e'^'^^ e"*"'" , 1, . . . , l) 



+M'^diag(l,...,l,e*'^"+S...,e*'^^)| , 6 = +,- 



(6.29) 
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Now the assumption ( 6.24 ) entails 



lim ri%t) 

t— »oo 



ML ■ 


■ ML 


• 


• 


ML ■ 




• 


• 






1 • 


• 


M'm ■ 




• 


• 1 



(6.30) 



whence ( |6.25| ) is clear. The assumption ( |6.26| ) gives rise to the same limit matrix as 



in ( |6.3[1| ), except that one should add yMf ,^_^i,... ,Mff^_^^j to the (n + l)th column. 
A moment's thought then yields 



lim r]Lt) 

t~*oo 



+ 



M: 



71+1 



(6.31) 



and so (|6.27D results. 

Proof of Proposition 6.1. We only show 
\1n + L{xo + sot,t)U~'^' 



lim 



exp 



(6.32) 



the proof for t 
entails 



In + L(xo + sot,t)U\ 

—oo being similar. Let us first note that the ordering of the positions 



< 



-v{a^ ) < u(a+_^) 



< 



< via 



(6.33) 



For So > v{af) we therefore have Aj(so) < 0, i 
Lemma 6.2 with 



M± = L(xo,0)C/^\ 
and n = in (|6.24| ) to obtain limit 1 on the Ihs of ( |6.32 ). Since we have 
fW(xo,so)=0, so>v{af), 



,N, cf. (6.20). Thus we can invoke 

(6.34) 



(6.35) 



we obtain ( |6.32| ) for sq > v{ai). 

For So S {v(a2),v{a'f)) we have Ai(so) > 0, with all other Ai(so) still negative. Thus 
Lemma 6.2 applies with (6.34) in effect, and with n = 1 in ( 6.241 ). Then the Ihs of ( |6.32 ) 
yields ex.p{2iaf), which equals the rhs (recall ( |6.18[) ). Likewise, it follows that ( |6.32D holds 
true for so not equal to the velocities ( |6.33| ). 

Next, we choose so = v{af). Then we have Ai(so) = 0, with all other Aj(so) negative. 
Thus we can exploit Lemma 6.2 with n = in ( |6.26| ), which yields a limit 



1 + (L(xo, 0)[/-i)ii _ 1 + exp{ef - 2afxo)V{^ exp{iaf 



l + {L{xo,0)U)u l + exp( 



2a]^xo)VL exp{—ia~[ 



(6.36) 



This agrees with the rhs of 06:321 ), cf. (0), and ( |6T0|). 

Choosing now so = f (a^). Lemma 6.2 with n = 1 in (|6.26| ) yields the limit 



exp(2ia]^ 



1 + 1712 exp(ia^ 
1 + m2 exp(— 



(6.37) 
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where m2 denotes the quotient of the principal minor of L(xq, 0) with respect to indices 1, 2 
and its 11-element. Cauchy's identity yields 

ma = exp(6'J - 2a^xo)F2^ tanh^[(g+ - q^)/2], (6.38) 

and so the limit agrees with the rhs, viz., 

exp [2iat -if+iqt,0t;xo + Atiq+,q-)/Aa+)] , (6.39) 



cf. ( |6lO| ), and ( pl^ ). 



Proceeding in the same way for sq equal to the remaining velocities, it is readily checked 
that ( |6.32 ) holds true as well. Therefore, Proposition 6.1 now follows. I 



Finally, we briefly consider the above A^-soliton solutions from the perspective of Sec- 
tion 7 in Ref. ||5|, especially as concerns the issue of soliton space-time trajectories. Inspec- 
ting the derivatives ( |2.55| ) and ( 2.6C| ), we see that there is an obvious choice of space-time 



trajectories for the 1-soliton solutions, namely, 

x+{t) = 6'+/2a+ + ?;(a+)t, (6.40) 

x-{t) = e-/2a~ - v{a-)t. (6.41) 

Let us now choose = N, = until further notice. Then we are dealing with the 
self-dual llrdiT = 7r/2) system. For the soliton solutions arising in Section 7 of Ref. 
this specialization amounts to the pure soliton case studied in Subsection 7A. In order to 
compare the present setting to loc. cit., we should take /3 = l,r = 7r/2in loc. cit., and 
substitute q, 6 ^ 9, q. Then the above Lax matrix (|5.2[) turns into a diagonal similarity 
transform of the Lax matrix (7.4) in loc. cit. 

Substituting next in loc. cit. 

y ^ X, aj = 2a^, Vj = v(a^), (6.42) 

one gets agreement with Eq. (7.7). Since L{x,t) ( |6^ ) amounts to a diagonal similarity 
transform of the matrix A{t, x) (7.6), we may apply Theorem 7.1 to the case at hand. Thus 
we obtain non-intersecting soliton space-time trajectories with long-time asymptotics 

x+ . (t) = (9+ T l^tiQ^)) + vi»t)t + 0(exp(Ttr,)), t ^ ±00, (6.43) 
where 

A+((Z+) ^ I E - E I 1^ (coth2[(g+ - g+)/2]) , qquadj = 1, . . . , iV, (6.44) 
\k<j k>jj 

rj = min2a'l\v{a'l) - v{a~![)\, j = 1,. . . ,N. (6.45) 

(See Ref. [^] for a picture of sine-Gordon soliton space-time trajectories.) 

To be sure, without prior knowledge it would not at all be obvious that the termi- 
nology used here is appropriate for the solution (|6.4|). Indeed, since L{x,t) is multiplied 
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by the g^-dependent phase matrices U and U* , we have no analog of the "1-particle su- 
perposition" formulas Eqs. (7.10)~(7.14) in loc. cit.. More importantly, at face value the 
long-time asymptotics of the spectrum of L{x, t) seems to have no bearing on the long-time 
asymptotics of the solution ^{x,t) (|6^). 

Even so, we need only invoke Proposition 6.1 for the case = N, N- = at hand 



to deduce that the asymptotic space-time trajectories in (|6.43D coincide with the 1-soliton 



trajectories in ( |6.9[ ), cf. also ( 6.40 ). Accordingly, the trajectories in ( |6.43| ) do encode the 



physical characteristics of the right-moving A'^-soliton solutions. 

Likewise, for the case A'^-i. = 0, N- = N we can use Theorem 7.1 in loc. cit. to obtain 
non- intersecting space-time trajectories that coincide with the obvious "locations" of the 
N left-moving solitons for asymptotic times. 

Whenever N^N^ > 0, however. Theorem 7.1 can no longer be used to define space- 
time trajectories. This is because L(x, t) need not have positive and simple spectrum when 
particles and antiparticles are present. Indeed, as (g"*", q~ ,6^, 9^) varies over Q ( 5.3S| ), the 



eigenvalues of L(0, 0) range over all of the right half plane, cf. loc. cit., Subsection 2C. 
As an illuminating example, consider the case A^-f, = 1 with q = (q~^ — q~)/2: 

T('rf\=( ta'^^l?! Vcoshg \ / exp(6i+(a;,t)) \ , . 

^'^ \i/coshq twh\q\ )'\ exp(^-(x, t)) ^ ' ^ ' 

e\x,t) = e^ -2a^ (^x-6v{a^)ty 6 = +,-. (6.47) 
Setting 

a{x,t) = tanh2([6l+(x,i) -e-{x,t)]/2) -l/cosh^g, (6.48) 

it is routine to verify that for a+ ^ a" the matrix L(x, t) has two distinct positive 
eigenvalues when (t(x, t) > and a complex-conjugate pair of eigenvalues in the right half 
plane when a{x, t) < 0, whereas L{x, t) is not diagonalizable when ^{x, t) vanishes. For an 
arbitrary fixed tQ, one therefore finds that L{x, to) has two distinct positive eigenvalues for 
X > Xq and x < Xq , and non-real spectrum for x G {xq ,Xq). Moreover, in the non-generic 
case = a~ the spectral character only depends on to. (As a bonus, these observations 
show that the Hamiltonians on 0, that generate the x- and t-flows generically do not leave 
the spectral decomposition of invariant.) 

In spite of this quite different state of affairs for A_|_A_ > 0, we believe that a more 
refined spectral analysis of L(x, t) should yield N non-intersecting space-time trajectories 
for t > T(+) and t < T^^\ where T^^^ depend only on {q, 6) G Q ( p. 39 ). The asymptotics 
of these trajectories should read 

x+ . {t)^-L.fe+TlA+{q+,q~))+v{apt, i ^ ±oo, (6.49) 
N+-j+l -^"i V ^ / 

X- ^ (t) ~ (ei T Wiq^, q-)) - v{ai)t, t ^ ±oo, (6.50) 



(with and A^ given by (6.10) and ( |6.11D ), so that the trajectories coincide asymptot- 



ically with the 1-soliton trajectories following from (6.9). 
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To explain what is involved here, it is illuminating to reconsider the example ( |6.46| )- 
( 6.48| ) with = —q~ (so that a"*" = a~). Then one readily verifies the following fea- 



tures. First, there exist unique T^+\ J'(-) such that L{x,t) has non-real spectrum for 
t G (T^~\T^^^) and distinct positive eigenvalues for t > T^'^^ and t < T^~\ Second, 
for t > T^"*") there exist unique x^{t), x~{t) with x~^{t) > x^{t) such that L(x^{t),t), 
6 = +,— , has non-degenerate eigenvalue 1. Likewise, for t < T^~^ there exist unique 
x'^{t), x~{t) with x'^{t) < x~{t) such that L{x^{t),t), 5 = +,—, has non-degenerate eigen- 
value 1. Finally, the asymptotics of the trajectories x^{t), 6 = +,— , is given by ( |6.49D 
and ( 6.50| ). (Physically speaking, the particle and antiparticle associated with this special 



two-soliton solution form a virtual bound state/resonance for t S [T^"-*, T*- "*"•*].) 

Of course, this example may be viewed as too special to yield convincing evidence. 
Indeed, a far more important reason why it is plausible that trajectories with the above- 
mentioned properties exist is the following lemma. 



Lemma 6.3. Let M e M and let 

E{t) = Me*^, (6.51) 



where D is given by ^6. 25 ) and ( 6.2(\ ). Then there exists T G M such that for all t > T 



the matrix E{t) has a non- degenerate eigenvalue en{t) obeying 

e„(t) = m„(l + p„(t)), mn = \Mn\/Mn-i\, (6.52) 

\Pn{t)\ < C exp{-trn), r„ = min((ii, . . . -dn+2, • • • , -^at). (6.53) 

Proof. This lemma can be obtained as a corollary of the proof of Theorem A2 in Ref. |21]. 
The latter theorem concerns the case where di, . . . ,dN are distinct, which is not implied 
by our assumption (|6.26| ). But when one follows the arguments in its proof, one easily 
sees that they entail the assertion of the lemma. I 

We now explain the bearing of this lemma on our trajectory conjecture. To this end 
we consider (for example) the matrix 

L{xo + v{a+)t,t) = L{xo,0)diag(exp{Xi{v{a+))t),...,exp{XN{v{a+))t^ , (6.54) 

where we used (H^) and ( |63o| ). Taking M = L(xo,0) and n = J in Lemma 6.3, we see 
that the assumptions are satisfied. Since the linear functions Aj(s) ( |6.20 ) intersect for s 



varying over a finite set I, the numbers di, . . . ,dN in the lemma are not distinct whenever 
v{aj) G I. This is the reason why Theorem A2 in Ref. ||2^ does not apply in general. 
(For s-values not in Z, however. Theorem A2 does apply to L{xq + st,t), entailing simple 
and positive spectrum for |t| sufficiently large.) 

As a consequence, we deduce that for t > Tj the matrix ( |6.54 ) has a non-degenerate 



eigenvalue ej{xo,t). The trajectory Xj{t) should now be defined so that 

ej{xj{t) -v{aj)t,t) = 1. (6.55) 

One readily checks that this would yield the desired asymptotics for t — > oo, but we cannot 
rigorously prove that there exists a unique xq € M such that ej{xo,t) = 1. (For one thing, 
the choice of Tj depends on the xq we fix.) 
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In any event, by now it will be clear what spectral feature our trajectory conjecture 
amounts to: For large times t, the matrix L{x, t) should have a non-degenerate eigenvalue 1 
for N and only N positions x~^_{t) < ■■■ < x^{t) < x^^{t) < ■■■ < xf{t) — the 
soliton space-time trajectories. Though the above arguments do not constitute a complete 
proof, they provide considerable evidence. To conclude, it should be emphasized that 
for N^N^ > one should not expect global trajectories. Indeed, as we have already 
established explicitly for the above example ( |6.46 ) with = —q^, it is likely that when 
left- and right-moving trajectories collide, the corresponding eigenvalue pair becomes non- 
real. 
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